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ABSTRACT 

We consider a general method of deprojecting two-dimensional images to reconstruct 
the three dimensional structure of the projected object, assuming axial symmetry. The 
method consists of the application of the Fourier Slice Theorem to the general case where 
the axis of symmetry is not necessarily perpendicular to the line of sight, and is based on 
an extrapolation of the image Fourier transform into the so-called cone of ignorance. The 
method is specifically designed for the deprojection of X-ray, Sunyaev-Zeldovich (SZ) and 
gravitational lensing maps of rich clusters of galaxies. For known values of the Hubble 
constant H and inclination angle, the quality of the projection depends on how exact 
is the extrapolation in the cone of ignorance. In the case where the axis of symmetry 
is perpendicular to the line of sight and the image is noise-free, the deprojection is 
exact. Given an assumed value of Hq, the inclination angle can be found by matching 
the deprojected structure out of two different images of a given cluster, e.g., SZ and 
X-ray maps. However, this solution is degenerate with respect to its dependence on the 
assumed H , and a third independent image of the given cluster is needed to determine 
Hq as well. The application of the deprojection algorithm to upcoming SZ, X-ray and 
weak lensing projected mass images of clusters will serve to determine the structure of 
rich clusters, the value of H , and place constraints on the physics of the intra-cluster 
gas and its relation to the total mass distribution. The method is demonstrated using 
a simple analytic model for cluster dark matter and gas distributions, and is shown to 
provide a stable and unique reconstruction of the cluster 3D structure. 

Subject headings: cosmology: observations - dark matter - gravitational lensing - dis- 
tance scale - galaxies: clusters: general 
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1. Introduction 

Uncovering the three-dimensional (3D) structure 
from observations of projected quantities is a funda- 
mental problem of astrophysics research. The 3D 
structure of rich galaxy clusters is particularly in- 
teresting, as the cluster geometry impacts on deter- 
minations of the Hubble constant from X-ray and 
Sunyaev-Zeldovich (SZ) measurements, cluster mass 
and baryon fraction determinations, and the underly- 
ing galaxy orbit structure. 

The most straightforward approach when dealing 
with projected quantities is to assume a spherical 
symmetry. The assumption is often applied within 
the framework of parametric models, for example the 
standard /3-model (e.g., Cavaliere & Fusco-Femiano 
1976; see Rephaeli 1995 for a review), or a non- 



parametric geometric approach (Fabian et al. 1981). 
Such simplifications can lead to biases in whatever 
quantity is being extracted from the data, especially 
when the systems under scrutiny might be biased to a 
certain geometry such as, for example, being detected 
since they are preferably elongated along the line of 
sight (hereafter LOS). This bias can be considerable 
as was illustrated in the particular case of merging 



clusters (Ftoettinger, Stone & Mushotzky 1997). 

Here, we restrict ourselves to the problem of re- 
constructing the 3D structure from maps of the X- 
ray emission or the SZ decrement (say), assuming the 
underlying 3D distribution is axially s ymmetric. W e 
show that the Fourier Slice Theorem ( Rybicki 1987 ), 
can be extended to general case of an arbitrary an- 
gle and the deprojection can be done by applying a 
simple approximation. In the case of astronomical im- 
ages the inclination angle is unknown. However, with 
the availability of SZ, X-ray and lensing maps, we 
can uniquely determine the cluster inclination angle 
by comparing the deprojections of the various images. 
This result enables us to ascertain the Hq dependence 
in the more general case where an axial symmetry is 
assumed. We apply this technique for the reconstruc- 
tion of the 3D cluster structure to simulated SZ and 
X-ray images in an analytical cluster model. 



2. Deprojection 

The deprojection of axially symmetric quantities is 



a classical problem (e.g., Lucy 1974), considered for 



determining the 3D stellar orbit structure in elliptical 



galaxies (e.g., Dchnen & Gerhard 1993; Dchncn & 
Gerhard 199^; [Binncy, Davies fc Illingworth 1990|), 



and in determining cluster 3D structure from X-ray 
observations (Fabricant et al. 1984). 

We deproject cluster images as follows. We adopt 
the convention that bold-face symbols denote 3D 
quantities (e.g., k — (k x ,k y ,k z )). Let the observer's 
coordinate system be defined with the Cartesian axes 
(x' , y', z'), with the z' axis aligned with the LOS. De- 
note the (cluster) source function coordinate system 
by the axes (x, y, z) where the z-axis is the cluster 
symmetry axis, forming an angle 9i with respect to 
the LOS. Let I(x',y') = J\'(x')dz' denote a pro- 
jected quantity (image) of the source function A. The 
3D Fourier transform (FT) of the source function is 
related to the image transform by 



e l-i(k' x *'+Ky')h(x',y')dx'dy' 



Ik' (k x , k y ). 



(1) 



The relations are easily seen in the case Qi = 90°: 
since the cluster has a rotational symmetry, in the 
cluster frame A(r) = A(r, z) and its FT is Ak(k) = 

A k (fc, k z ) where k = ^jk x + k%. Now k = ^Jk'? + k'y 2 

and k z = k' x , and the relation between the FT of the 
source function in the observer and cluster rest frames 



is A k , (k') = Ak yyky + k' 2 , k' x J . Deprojection is ob- 
tained by inverse transforming, and using equation 
(III), to obtain 



A(r,z) 
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(27T)2 



e lk ' z I k ,(k,k z )J (kr)kdkdk z . (2) 



This particular case where the symmetry axis is 
perpendicular to the LOS is an exact application of 
the Fourier Slice Theorem. The general case is easy to 
obtain. For arbitrary inclination angles, equation (Q) 
holds. The expression relating the FT of the source 
function in the observer and cluster rest frames is ob- 
tained simple by coordinate rotation, so that 

A' fc ,(k') = A fc (y(-fci sm9 t + k' x cosfl,;) 2 + kg, 



k' z cos 9i + k' x sin 9, 
\ k (k 7 k z ). 



(3) 



where the last equality again holds due to axial sym- 
metry. Inverse transforming, we find the desired ex- 
pression for the source function in real space 
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A doproj (r, z) 



(2tt) 2 



kdkdk z e tkzZ J (kr) x 



sin ( 



k 2 — k 2 cot 2 I 



(4) 



For wave- vectors k < \k z \ cot Qi the argument of the 
image FT becomes imaginary. This defines what has 
been dubbed the "cone of ignorance" (iRybicki 198?| ; 
Gerhard & Binney 199C ) inside which there is no in- 
formation on the 3D structure in the image. Gener- 
ally, A depr °j(r, z) 7^ A(r, z), and the difference between 
the two depends on the power hidden in this region. 
For clusters aligned along the LOS, this region com- 
pletely covers k-space, and the method fails. In the 
other extreme, Qi = 90°, all of the information re- 
quired for a unique deprojection is contained within 
the image. 

The implementation of this method is straight- 
forward, assuming for the moment that the inclina- 
tion angle is known. The image FT is evaluated at 
wavevectors 



k' x = k z / sin Qi and k' y = ^ k 2 - fc 2 cot Q 2 (5) 

and the inverse transform applied. Numerically, this 
requires some interpolation between the grid points of 
a finite FT. For moderately smooth underlying clus- 
ter distributions, this should not pose difficulties al- 
though this can be rigorously tested with numerical 
cluster simulations, for example. We consider a sim- 
ple application in §^|. 

The cone of ignorance is a problem for the unique- 
ness of the deprojection. Palmer (1994) shows that 
any 3D axisymmetric source function can be uniquely 
recovered via deprojection if it can be expanded as a 
finite sum of spherical harmonics. However, Gerhard 
& Binney (1996) have constructed an example that 
obeys this theorem and lurk only within the cone of 
ignorance, while still having a physically interesting 
density distribution (a disk like structure). This type 
of degeneracy has limited the utility of this technique 
in recovering the 3D density distribution in galaxies, 
where this type of structure is plausible. However, 
Romanowsky & Kochanek (1997) employed a numer- 
ical regularization method deal with information hid- 
den in the cone of ignorance, and constrained physi- 
cally acceptable models by comparing with observed 
stellar velocity dispersion. 



Clusters are more favorable systems for this method. 
Often cluster images are elongated along an (approx- 
imate) axis of symmetry. In a bottom-up clustering 
scenario, cluster formation proceeds by an ongoing 
merging process, and one expects clusters not to be 
highly flattened, and to be preferentially elongated 
objects. Moreover, clusters form via dissipationless 
clustering and therefore lack sharp discontinuities and 
irregularities, resulting in a relatively smooth and reg- 
ular Fourier space behavior. Given this, at least some 
of the power resides outside the cone of ignorance, 
which permits extrapolation into this region. More- 
over observations of the weak lensing effect, the X-ray 
surface brightness and the SZ decrement provide pro- 
jections of different powers of the gas and mass den- 
sity distributions so that only very contrived models 
will be degenerate in all three images. The idea is 
then that by combining X-ray, SZ and lensing ob- 
servations, we can hope to determine the inclination 
angle by comparing the deprojection of the various 
images - a variant of this idea was considered some 
time ago in an analysis of the X-ray emission and the 
galaxy distribution in the cluster A2256 (Fabricant e\ 



1984) 



3. Noise properties 

To compare the deprojection from various wave- 
lengths statistically, we need to characterize the noise 
properties of the source functions. As with all de- 
projection methods, the noise is correlated in the re- 
sulting 3D map. Suppose the image is pure white 
noise, with rms a, and hence (I(x'i, y[)I(x , 2, y 2 )) = 
(2ir) 2 5(x[ — x 2 )S(y' 1 — y' 2 ) a 2 . The noise correlation 
function in the deprojected image is given by 



(A(ri,zi)A(r 2 ,z 2 )} = 0' 



sin 9i 



kdKdke lk ^ Zl ~ z ^ x 



y/k 2 - hi cot t JoQsri) J (fcr 2 ). (6) 

Not surprisingly, the high frequency components dom- 
inate and hence some smoothing must be employed. 
The analytic nature of the correlation matrix is an 
attractive feature of the method. This allows a direct 
statistical comparison of the source functions deter- 
mined independently from SZ and X-ray observations, 
and alleviates the need for simulations to determine 
the statistical significance of features found in the de- 
projections. 

A further nice feature of the method is that any 
isotropic smoothing function applied to the image will 



3 



depend only on the modulus of the radial wavevec- 
tor, |k|, in the deprojection. This can easily be seen 
since the identification of coordinates between the ob- 
server's frame and the frame of reference of the source 
function given in equation (||) yields k' 2 — > k 2 + k 2 = 
k 2 |. A simple and practical Gaussian smoothing fil- 
ter gives a particularly simple correlation of the noise 
(A(r) 2 ) txr" 1 . 

4. Practical Application 

To illustrate the method, we consider a simple an- 
alytical model for the cluster gas distribution. We use 
a simple gas density profile, p g (x,y,z) oc [r 2 + (x 2 + 
y 2 ) I a 2 + z 2 1 c 2 ] -1 and create simulations of SZ and X- 
ray observations for a prototype cluster with the fol- 
lowing characteristics: r c — 100/i _1 kpc, T — 7 keV, 
central density of p g o — 0.01 cm~ 3 , axis ratio of 3:1 
and an inclination angle of 0, = 60° . White noise is 
added to mimic levels of ROSAT PSPC observations, 



and SZ maps of Carlstrom, Joy & Grego (1997) 



The procedure then to determine the cluster 3D 
structure is as follows: given X-ray and SZ images, we 
apply a smoothing filter to reduce the noise levels in 
the images. Motivated by the high degree of noise re- 
moval necessary for this type of application, one of us 
(Hoffman 1998) has proposed a method for noise re- 
moval by wavelet thresholding. This method has the 
advantage of superior noise suppression over simple 
linear filtering, while maintaining spatial resolution. 
After denoising, an inclination angle is assumed, and 
the FT of the images formed. To perform the depro- 
jection, we employ equation (Q), and evaluate the FT 
at wavevectors given by equation (||). This requires 
some interpolation from the finite image FT, and we 
employed cubic splines. 

For smooth cluster density profiles, the image FT 
is quite regular and the extrapolation into the cone of 
ignorance can be done robustly. For this example, we 
employ a relatively simple linear extrapolation, with 
the slope being taken from the tangent to the con- 
tours at the cone boundary. If significant substruc- 
ture exists, the worry is that this sort of extrapolation 
would smooth over the structure. Clearly, more so- 
phisticated approaches could be adopted as the data 
warrants, however this simple extrapolation seems to 
work well for relatively smooth cluster profiles. We 
are currently undertaking a more thorough investi- 
gation using cluster SPH simulations, and we expect 
that the simplistic approach adopted here will not be 



the final word on the subject. 

Having evaluated the image FT at the wavevectors 
in equation (j^) and extrapolated into the cone of igno- 
rance, we are in position to perform the deprojection. 
If the inclination angle is known, the reconstructed 
3D structure of the source functions \sz(r, z) and 
Xx(r,z) is excellent. Of course the fact we have ig- 
nored heretofore is that the method requires a priori 
that the inclination angle is known, whereas this is 
not the case for astronomical observations of clusters. 
For any one image, all one can do is assume some 
angle. 

The error made assuming the incorrect inclination 
angle is displayed in Figure [l]. The top left and the 
bottom left and right panels show the deprojection 
of the SZ image compared with the actual SZ source 
function for three assumed inclination angles. In the 
case where g^ ssumcd jL Q it a substantial discrepancy 
between the actual and deprojected source function 
is found. This is to be compared with the good 
agreement found when the correct 0i is assumed. To 
further investigate the robustness of the method the 
following comparison was made (Figure |l|, top-right 
panel): using the correct inclination angle, the depro- 
jected SZ source function had been used to to calcu- 
late the X-ray source function, which was then com- 
pared with the actual one. Again, very good agree- 
ment is found. 

This comparison of the various observables is ad- 
mittedly simplistic and the procedure for real obser- 
vations will necessarily be more involved. For exam- 
ple, the current best resolution and highest signal-to- 
noise SZ measurements being obtained at BIMA and 
OVRO are interferometric (Carlstrom, Joy & Grego 
1997). Since the SZ data is in Fourier space, with 
irregular spacing, the algorithm as presented here of 
deprojecting two or more images does not strictly ap- 
ply. However, it is easy to compare the X-ray and 
SZ data in some other fashion by, for example depro- 
jecting the X-ray image, constructing the \sz(^x) 
source function from the X-ray deprojection, project 
this source function into the SZ observational plane, 
and there make the comparison. This type of applica- 
tion will need to be tailored to the particular data set 
in hand and tested by numerical cluster simulations. 

Note that weak lensing observations provide a fur- 
ther probe into the cluster 3D structure. However, 
to compare the observable from weak lensing, the to- 
tal mass surface density, with X-ray and SZ observa- 
tions, which probe the gas density and temperature 
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Fig. 1. — Reconstruction of A(r) at fixed z for dif- 
ferent inclinations, showing actual and reconstructed 
images. The top left, and bottom left and right panels 
show the SZ deprojection (solid lines), compared with 
the true underlying source function (dashed line), as- 
suming three different inclination angles. The top 
right panel shows the gas density inferred from the 
X-ray and SZ maps independently, with the correct 
value of the inclination angle used in both cases. The 
agreement in the X-ray and SZ deprojections is quite 
good, if the inclination angle is known. This suggests 
this is a viable method for determining the cluster 
geometry by requiring consistency between the two 
deprojections. 



distributions, we need to appeal to an assumption for 
the cluster dynamical state. Assuming the gas is in 
hydrostatic equilibrium in the cluster potential well, 
the gas pressure, P, is related to the potential, 0, by 
VP = —p g V(f)(r) and the total mass density is deter- 
mined via the Poisson equation V 2 (f> = AirGp. This 
illustrates the procedure for comparing the maps: 
given SZ and X-ray images, one deprojects to deter- 
mine the gas density. From weak lensing observa- 
tions, one determines the total mass density. Compar- 
ing the deprojections by assuming hydrostatic equilib- 
rium will then yield further insight into the dynamical 
state of the cluster, and the relative distributions of 
the dark matter and the gas. 

5. Cosmological Parameter Dependence 

The deprojection of images of objects at cosmolog- 
ical distances is further complicated by the Hq depen- 
dence, but it also provides an interesting and practical 
way for determining the value of Hq. We first recall 
the Hubble constant dependence on the X-ray and SZ 
maps. Let the observer's coordinate systems be de- 
noted by (9 x ,9 y ,z), where z is the cluster redshift. 
Physical coordinates are given by x! = Da(z)0 x , 
y' = Da{z)9 v where Da(z) is the angular diameter 
distance to the cluster. The Hubble constant depen- 
dence can be factored out as Da(z) oc h^ 1 and hence 
we rescale the variables by x! = h x' (similarly for 
y). Expressing the image and the object structure in 
these coordinates, the Hubble constant dependence is 
explicit: the 2D FT of the image yields an h~ 2 scaling 
and the 3D inverse FT gives an h 3 scaling. It follows 
that the expression for X(r, z) (equation |^) varies as 
h. 

We can see this explicitly, adopting a model for 
the intracluster gas. Let V to be the projection op- 
erator and P -1 be the inverse. The X-ray source 
function (which is proportional to p 2 g T x l 2 ) is given 
by X x (r,z) = hP^ 1 [I x (9 x ,9 y )] (and similarly for 
the SZ source function, X sz cx p g T). Then given 
a model that relates X x and A S2: , one can solve for 
h. For example, for a polytropic gas (p oc and 
T oc pP, P = 7 — 1), one finds 



i - ' 

2-/3 



h = 



(p-i[4])^ 



(7) 



In the model employed here, the deprojection oper- 
ator scales with sin^,;. Thus V7 1 oc sin (9, , and there- 
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fore each of the source functions from the X-ray and 
SZ deprojection has a scaling with h sin 9i. This is 
illustrated in Figure |2|. We created simulated noise- 
free SZ and X-ray maps, with the density model de- 
scribed in §^ and a polytrope of 7 = 5/3. We applied 
the deprojection to both images. We took the depro- 
jected Xsz, predicted the X-ray source function, A^, 
and compared with the deprojection of the X-ray im- 
age. In the three panels of Figure || we display the 
results for three assumed inclination angles. In all 
cases, the shape agreement is almost perfect, but the 
normalization disagrees. 

The relation between the true and assumed incli- 
nation angle and the true and inferred Hubble con- 
stant is shown in the top right panel of Figure ||. The 
basic result is that the value for Ho inferred scales 
with sin6>i. We have not succeeded in proving this 
h sin 9i degeneracy in the general case of an arbitrary 
3D structure, but we speculate that this behavior 
characterizes (at least to first order) cluster images 
in general. 

The hshiOi degeneracy can be interpreted either 
as a positive or negative feature of the method. The 
encouraging fact is that, at least in the context of the 
model of the cluster gas used here, requiring agree- 
ment between the deprojected source functions will 
place joint constraints on Hq and the inclination an- 
gle. The disappointing aspect is that from these two 
maps alone, one can not constrain the cluster geome- 
try and Hq independently. 

6. Conclusion 

We present a new algorithm for constructing the 
3D structure of axisymmetric objects from projected 
images. The method is shown to be especially suited 
to the deprojection of images of clusters of galaxies, 
and is deonstrated to be stable in the presence of noise 
and quite a 'wide' cone of ignorance. It is shown that 
given an assumed value of Hq and a model relating 
the gas density and temperature, the cluster inclina- 
tion angle can be solved given its SZ and X-ray maps. 
With another image of the cluster, Hq and the incli- 
nation angle can be determined independently. The 
application of the method to realistic cluster simula- 
tions, and to SZ, X-ray and weak lensing observations, 
will be presented in forthcoming papers. 
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Fig. 2. — The effect of the (unknown) inclination an- 
gle on the determinations of Hq. The bottom and 
top left panels compare the X-ray source function 
(oc p 2 T 1 / 2 ) in dashed contours, with the prediction 
from the deprojection of the SZ map, using 3 differ- 
ent assumed inclination angles. Good agreement is 
obtained in the shape of the X-ray isophotes for any 
assumed inclination angle. However, the amplitude is 
mismatched and this translates in an error on the in- 
ferred value of -ffo, shown in the top right panel. The 
inferred value of Hq scales with the inclination angle 
as sin#i. 
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